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On Convergence Rate of Leader-Following Consensus of 
Linear Multi-Agent Systems with Communication Noises 

Long Cheng, Yunpeng Wang, Wei Ren, Zeng-Guang Hon, Min Tan 


Abstract —This note further studies the previously proposed eonsensus 
protocol for linear multi-agent systems with communication noises in 1151 . 
(S). Each agent is allowed to have its own time-varying gain to attenuate 
the effect of communication noises. Therefore, the common assumption 
in most references that all agents have the same noise-attenuation gain is 
not necessary. It has been proved that if all noise-attenuation gains are 
infinitesimal of the same order, then the mean square leader-following 
consensus can be reached. Furthermore, the convergence rate of the 
multi-agent system has been investigated. If the noise-attenuation gains 
belong to a class of functions which are bounded above and below 
by {0 G (0,1)) asymptotically, then the states of all follower 
agents are convergent in mean square to the leader’s state with the 
rate characterized by a function bounded above by asymptotically. 

Index Terms —Multi-agent systems, leader-following consensus, noises; 
time-varying gain; convergence rate. 

I. Introduction 

Communication noise is an unavoidable factor in the distributed 
consensus of networked multi-agent systems. It has been found in 
m that the traditional protocol cannot solve the consensus problem 
with the existence of the communication noise. Therefore, how to 
effectively attenuate the noise effect becomes an interesting research 
topic. One popular idea is to employ a time-descending gain (some¬ 
times called the stochastic-approximation type gain) in the consensus 
protocol. In the early study phase, many scholars investigated the 
stochastic consensus of first-order integral multi-agent systems with 
communication noises. For example, under the fixed topology case: 
the stochastic-approximation type gain was first employed to solve 
the mean square and almost sure consensus problems (3; Li and 
Zhang proved that the mean square consensus can be achieved in 
the continuous-time domain if and only if the noise-attenuation gain 
satisfies the stochastic-approximation type condition jT); some nec¬ 
essary and sufficient conditions for ensuring the stochastic consensus 
with both communication noises and time delays in the discrete¬ 
time domain were presented in (3. Under the switching topology 
case, the stochastic-approximation type protocols were also proved 
to he effective in both the discrete-time domain a-E] and the 
continuous-time domain (6). It is also noted that for the continuous¬ 
time mean square leader-following consensus problem, the necessary 
and sufficient conditions of noise-attenuation gains can be slightly 
relaxed compared to the leaderless case Q, (U The aforementioned 
papers study the additive noise while the multiplicative noise has 
also been considered in ©-[m It is usually assumed that the noise 
intensity is proportional to the state differences between agents. In 
particular, if the state difference becomes zero, then the noise effect 
disappears. Therefore, the protocol for solving multiplicative noises 
may not require the time-descending gain. A few recent results also 
provided the protocols for higher-order integral multi-agent systems 

Pai't of this work has been presented at the IFAC World Congress, 
Capetown, South Africa, August 2014. L. Cheng, Y. Wang, Z.-G. Hou and 
M. Tan are with the State Key laboratory of Management and Control for 
Complex Systems, Institute of Automation, Chinese Academy of Sciences, 
Beijing 100190, China. W. Ren is with Department of Electrical and Computer 
Engineering, University of California, Riverside, CA 92521, USA. 

Please address all correspondences to Dr. Long Cheng at Email: cheng- 
long@compsys.ia.ac.cn; Tel: 8610-82544741; Fax: 8610-82544794. 


with communication noises. For instance, it has been proved in fl2l 
that the stochastic-approximation type gain is still the necessary and 
sufficient condition for ensuring the mean square average consensus 
of second-order integral multi-agent systems under the fixed topology. 
Further results regarding the switching topologies were presented 
in HU, (H. For the general linear multi-agent system, there are 
also some attempts on solving the mean square consensus and 
the almost sure consensus under the fixed topology 03 and the 
switching topology dj. However, there are still some limitations 
in the current study of higher-order integral multi-agent systems 
with communication noises. For example, all agents are required to 
have the same noise-attenuation gain in most existing publications, 
which implies that certain global information should be known by all 
agents; and the convergence rate of the multi-agent system is rarely 
considered. 

Motivated by the above observation, this note first modifies the 
consensus protocol for general continuous-time linear multi-agent 
systems proposed in dsi. In the modified protocol, each agent is 
equipped with its own noise-attenuation gain. Hence there is no 
need to assume that all agents have the same gain. It is proved that 
the mean square leader-following consensus can be reached by the 
modified protocol if all noise-attenuation gains are infinitesimal of 
the same order. It is interesting to find that for the leader-following 
consensus of general linear multi-agent systems, the stochastic- 
approximation type requirement on the noise-attenuation gains can 
be relaxed (the square integrable condition is not necessary). Next, 
the convergence rate of the multi-agent system under the modified 
protocol is presented. Although the state-transition matrix has been 
explicitly obtained in G5l (under the assumption that all agents 
have the same noise-attenuation gain), the entire dynamical behavior 
of the multi-agent system can be determined by calculating the 
solution to the governing stochastic differential equation. It is still 
difficult to clearly tell the convergence rate of the multi-agent 
system since the solution to the governing stochastic differential 
equation is a very complicated function of noise-attenuation gains 
and the graph Laplacian matrix. Therefore it is more challenging 
to answer what is the convergence rate of the multi-agent system 
with the consideration of agent-dependent gains. Fortunately, by some 
recent results in Oil: if we assume that the noise-attenuation gain 
belongs to certain representative class of functions, then quantitatively 
determining the convergence rate becomes possible. In this note, if 
the noise-attenuation gain belongs to a class of functions bounded 
above and below by (0 G (0,1)) asymptotically, then the 
states of all follower agents are convergent in mean square to the 
leader’s state with the rate characterized by a function bounded 
above by asymptotically. This convergence rate analysis is the 
main improvement of this note compared to the previous conference 
version C3- 

Following notations are used throughout this note: C, R, N”*" 
denote the field of complex numbers, the field of real numbers 
and the set of positive integer numbers, respectively; In denotes 
the n-dimensional unit matrix; = (1, • • • ,1)^ £ On = 
(0, • • • ,0)^ G R"'; ® denotes the Kronecker product; for a given 
matrix X, denotes its transpose; ||A ||2 denotes its 2-norm; 


2 


for a random variable/vector x, E{x) denotes its mathematical 
expectation; For any two functions /(f) and g{t), f{t) = 0 {g{t)) 
represents limsup^^^ |/(f)/c?(f)| < oo; f {t) = o{g{t)) repre¬ 
sents limt->oo \f{t)/g{t)\ = 0; f{t) = 0 {g{t)) represents 0 < 
liminft^oo |/(f)/c/(f)| < limsup^^^\f{t)/g{t)\ < oo; for any 
X € C, Si(x) denotes its real part. 

11. Preliminaries 

Consider a multi-agent system composed of + 1 agents labeled 
from 0 to N. The communication network among agents is modeled 
by a digraph (J = {Vg, £g, Ag} where Vg = {vq, • • • , Vjv} denotes 
the node set; £g = {cij} C Vg x Vg denotes the edge set; 
and Ag = [dij] € is the weight matrix. Here Vi 

represents the ith agent; dj = {vj,Vi) € £g means that there is 
an available communication link from agent j to agent i\ aij > 0 
is the communication quality associated with the edge eij and 
dij > 0 if eij £ Eg, dij = 0 if Cij ^ Eg. Agent j is called 
the parent of agent i if dj £ Eg. The neighbor set of agent i is 
the set of all its parent agents, i.e., A/i = £ Eg}. The 

agent is called the leader if its neighbor set is empty, otherwise 
the agent is called the follower. The Laplacian matrix of Q is 
defined as £ = diag( aoi, • • • — Ag. There is 

a directed path from node nij to node Vi„ if there is a set of nodes 
{vi2 , • • • , } such that edges d^i ^, ■ • ■ , all belong to 

Eg. The digraph Q is called to have a spanning tree if there exists at 
least one node such that there are directed paths from this node to any 
other nodes in Q. In this note, it is assumed that the communication 
graph of the multi-agent system has a spanning tree. We assume that 
this multi-agent system has one leader (the leader is labeled by 0 ). In 
other words, the multi-agent system has a leader-following structure. 
It is obvious that the Laplacian matrix of such a multi-agent system 
has the following form 


agent i and its neighbors’ information such that limt_>oo E\\xi{t) - 
®o(f)|P = 0 and limsupj_,,j^ E\\xi{t)\\'^ < oo, Vi = 0 , 1 , • • • ,N. 

In this note, the leader-following consensus protocol for the ith 
agent is proposed as follows 

Ui{t) = KiXi{t) + ai{t)y^ aijK2{yij{t)-Xi{t)),i = 0 , - ■ ■ ,N, 

( 3 ) 

where di{t) > 0 is the consensus gain for the ith agent; Ki — 
(-01,-02 - for, • • • , -On - fon-i) and K2 = (fei, • • • , &n-l, 1 ); 
(61, • • • , bn-i) are parameters to be determined later. It is easy to 
see that the proposed protocol is different from the ones proposed in 
dsi, m because each agent has its own consensus gain di{t). 

III. Main Results 

Let Xpit) = (sf (f), • • • , xJi{t)}'^. Then substituting (O into l |3 
obtains that 

Xp{t) = {In ® (A -f BKi) - Ql{t)L2 ® BK2)XF{t) 

— 2l(f)Z/i ® BK2Xo{t) + 2l(f)Sr;(f), 

where 21 (f) = diag(ai(f), • • • ,dN{t)), E = diag(Ei,'-' ,Ejv), 
Ei = BK2{pio,pii,--- ,piN) and 77(f) is the nN{N + 1 )- 
dimensional standard white noise vector composed of rjij{t), i,j = 
OA,---_,N. 

Let Xpit) = Xpit) — Ijv 0 xo{t). Then 

Xpit) = (In (E> {A + BKi) -^{t)L20 BK2)XF{t) + ^{t)T.ri{t). 

Let X{t) = {In ®K2)XF{t). It is easy to see that K2{A-\-BKi) = 
0 % and K2BK2 = K2. Therefore, 

i(f) = - 2 I(f)L 2 X(f) -b {In ® A' 2 ) 2 l(f)E 77 (f). ( 4 ) 

Let $(f, to) denote the state transition matrix of H(f) = 
—2l(f)I/2H(f). By Ito integral formula, the solution to Q can be 
written as 


Lemma 1 (Theorem 2 in |[8)). If the communication graph Q of the 
multi-agent system has a spanning tree, then all eigenvalues of L2 
have positive real parts. Furthermore for any diagonal matrix D with 
positive diagonal elements, all eigenvalues of DL2 have positive real 
parts as well. 

The dynamics of the fth agent is described by 

Xi{f) = Axi{t)Bui{f), i = ( 2 ) 

where Xi{t) £ B"' is the state vector, Ui{t) G E is the control input. 
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X(f) = Ji(f,fo)-b J2(f, fo), ( 5 ) 

where Ji(f,fo) = 'l>(f, fo)X(fo), J2(f,fo) = J^^^{t, s){In 0 
A'2)2I(s)EdVL(s) and W{t) is the 7rA^(A-bl)-dimensional standard 
Brownian motion corresponding to 77(f). 

Throughout this note, the following four conditions hold. 

(Cl): d{f)dt = 00, where d(f) = maxi=i_... ,N{di{t)}. 

(C 2 ): There exist positive constants p.i < p.2 < 00, T < 00 and 
P € ( 0 , 1 ) such that for Vf > T, < ai{t) < p,2t~^, 

,N. 

(C 3 ): All consensus gains {ai(f), • • • ,ajv(f)} are infinitesimal of the 
same order as time goes to infinity. 

(C 4 ): All roots of the following polynomial have negative real parts 

S -b bn—iS -b • • • + b2S -b III = 0 . (6) 


(ai,--- ,a„) are coefficients determined by the agent’s essential 
dynamical characteristics. It is obvious that any controllable single 
input single output system can be transformed into this Luenberger 
canonical form. 

The agents exchange information via a noisy communication 
network. The information which the ith agent receives from the jth 
agent is denoted by yij{t) = Xj{i) -b Pijriij{t) where r]ij{t) = 
( 77 iji(f), • • • ,r)ijn{t))'^ £ E" is the 72 -dimensional standard white 
noise; pij = diag(piji, • • • ,pij„) G {\piji\ < 00, I = 

1, • • • , 72 ) denotes the noise intensity matrix. It is assumed that riiji (f) 
(f, j = 0,1, • • • , N; I = f, ■ ■ ■ ,21) are all mutually independent. 

The control objective is to achieve the mean square leader¬ 
following consensus. That is: design control inputs Ui{t) by using 


Since ai(f), • • • ,ajv(f) are infinitesimal of the same order as time 
goes to infinity, there must exist N positive constants ci,-- - , cjv 
such that limt_nx3 di{t)/d{f) = d, i = 1 , - ■ ■ ,N. 

Theorem 1 . If Conditions \(Cl)^C 4 ^ hold, then the proposed pro¬ 
tocol defined by 0) can solve the mean square leader-following 
consensus problem of 13 . Furthermore, the convergence rate of 
the multi-agent system is characterized by \\E{xi{t) — a;o(f))||2 = 

0{e ^ and E\\xi{t) — a:o(f )||2 = 0{t~^) {i = 

1, • • • , N), where pi > 0 and P £ ( 0 , 1) are defined in |(C 2 )[ 
Amin = min{5R(Ai), • • • , SR(Aiv)} > 0 , Ai, • • • , \n are eigenvalues 
of CL 2 and C = diag(ci,-- - ,cn); and e is any constant in 
( 0 , min{l, A min} )• 
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Proof: Since Q has a spanning tree, by Lemma □ Amin > 0 . 
It is easy to see that 2 t(i)I /2 = a{t)CL2 + a{t)D{t)L2, where 
D{t) = - ,djv(l)) and di{t) = ai{t)/a{f) — a. Then 

l^mt-^oo D{t) is a zero matrix. Therefore, by Lemmas and [ 7 ] for 
Ve G ( 0 , min{l, Amin}), there must exist two positive constants 
Ml, M2 < 00 such that for \ft > to. 


||$(f,fo)||2 < Mie ^*0 < M2e 1-? 

Therefore 


( 7 ) 

(8) 


Moreover, it is easy to see that E[J2(t,to)) = Ojv- Therefore, Vi = 
I,... ,iV, 

E{xi{t)) = E{K2{xi{t) - xo{t))) = 0 {e -‘ ), ( 9 ) 



Fig. 1. The trajectories of the state differences between the follower agents 
and the leader agent under the proposed protocol. 


which together with Lemma leads to the fact that \\E{xi{t) 

3 ;o(f ))||2 = 0 (e ^ f ^). 

It can be calculated that 

\\E{J2{t,to)J^{t,to))\\2 

= [ -L(Ls)(/iv®iT2)a(s)EE^a^(s)(7^® 

JtQ 

<iij;v®i^2ii^iiE|i^ r \mt,s)\\i\ms)\\Us 

Jto 

Jto 

which implies that 
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7?l|72(Lfo)||2 = (!7(||7?(J2(Lto)J2^(Lfo))||2) 


= O 






Fig. 2. The trajectories of the norms of all agents’ states under the leaderless 

( 10 ) 


By the same procedure of Lemma A.2 in HU, it can be proved 
that ^ 

Jto 

Hence, i7||X(f)||2 = which indicates that i?|7L2(a;i(f) — 

®o(f))|^ = i = 1, • • • , A^. By Lemma[^and Condi tion |(C4)| 

it can be obtained that i?||a;i(t) — a:o(f )||2 = i = 1, • • • ,N. 

Finally, by 0), the closed-loop dynamics of the leader agent is 
xo{t) = {A-\- BKi)xo{t). Since Condition |(C4)| holds, there must 
exist a constant v such that limt_>oo xo(t) = x* = (v,0, ■ ■ ■ , 0)^ € 
R". Therefore, 7?||a;i(f)||2 < 00 , i = 0, ■ ■ ■ , N. ■ 

Remark 1. Compared to the leaderless consensus studied in flTl . it is 
interesting to see that /3 in |(C2)| can belong to (0, 0.5), which means 
that the square integrable condition on ai{t) ai{s)ds < 00 ) is 
not necessary. From this point of view, the leader-following consensus 
seems easier to be achieved than the leaderless one. 


state differences between the follower agents and the leader agent 
are mostly upper-bounded by 5/t^'^, which is consistent with the 
theoretical analysis on the convergence rate. 

Next, a simulation example is conducted to show that the leaderless 
consensus studied in EH. (m needs the square integrable condition 
on the noise-attenuation gain. To this end, an extra edge eo4 is added 
(ao 4 associated with eo 4 is set to be 1), which results in a multi-agent 
system without any leader. In this case, p cannot be set to be any 
value belonging to (0, 0.5). For example, we set all noise-attenuation 
gains as the same value apt) = 1/(1 +f)^'*, i = 0, • • • , 4 . The 
simulation result is given in Fig. By the definition of leaderless 
consensus in 051,03, there must exist a random vector x* satisfying 
i?||a;*||^ < 00 such that limt_>oo i?||a;i(f) — a;* ||^ = 0, i = 0, • • • , 4 . 
This definition implies that limt_>oo E\\xi{t)\p < 00, i = 0, • • • , 4 . 
Although by Fig.|2 it seems that limsupj_,,^ i?||a;i(f) — Xj{t)\\^ < 
00 , the requirement that limt_>oo i?||a;i(t)||^ < 00 (i = 0, • • • , 4 ) 
cannot be satisfied. 


IV. Simulation Examples 


Consider a multi-agent system composed of five agents. In 
Xi{t) G R^, ai = —1, a2 = 1, 03 = 0 and 04 = —2. The elements 
of the weight matrix Ag are set as: aio = 020 = <124 = a4i = (143 = 
1, ffl 3 i = 2 and all other elements are zero. The noise intensities 
Piji in (ID are all set to be 1. The controller parameters in ([D are: 
Ki = (1, -2, -3, -1), K 2 = (1, 3, 3,1), ao{t) = 0.15/(f -f 1)° \ 
ai(f) = 1.2/(f -f 1)°-^, 02 ( 1 ) = 1.5/(3f -f 1)° '‘, 03 ( 1 ) = 0.6/(t + 
2)°'^ and 04(1) = 1.5/(4f-|-1)° '^. From the simulation results shown 
in Fig. [T] it can be seen that the leader-following consensus can be 
achieved in the mean square sense. In addition, the trajectories of the 


V. Conclusions 

This technical note relaxes the assumption made in 03 , lfi 3 that 
all agents have the same noise-attenuation gain. Each agent is allowed 
to have its own time-varying gain function. It is proved that if all 
consensus gains are infinitesimal of the same order, then the modified 
protocol can still solve the mean square leader-following consensus 
problem of general linear multi-agent systems. In addition, this note 
presents the convergence rate of the multi-agent system when the 
noise-attenuation gains belong to a representative class of functions. 
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Appendix 

Lemma 2. If Condition \(C2)\ holds, then V6 > 0, e = 

-bAtljl-/3 _o^ - 

0{e 1-^* ) = o(t where /3 € (0,1) A defined in Condition 

Proof: If Condition |(C2)| holds, then d(f) > Therefore, 


o *'Jtn ^ g ^'/tn PI® ^“*0 


= e“T^^ 


bMl . 


= C>(e"T^ 


w/iere A £ C, PQ{t,to) = e ant/ P/'(f,fo) = 

“ /t*o «('P)^i^-i(P to)Po{t, T)dT, i = 1, 2, ■ • ■ , r - 1. 


Lemma 5. For Ve > 0, A £ C, K(A) > 0 f/iere exAfj a positive 
constant M < 00 such that ||'l>A(f, fo )||2 < Me 
where is defined by ( 1131 . 


Proof: It is easy to see that |Po^(Lio)| = e .Aq 
A ssume that |P/'(t,fo)| < '^*0 where Mi is 

a finite positive constant. Then 


Let X = then Since /3 > 0, there must exist 

a positive integer n such that < n. By L’HospotaTs rule, it is 
obtained that 


lim 

t—^OO 


he-i jl-e b/J-i „ 3 

p 1 — 0 np 1 — 0 np 

= lim -H— = lim —r - < lim 


e~T^ 


t-P 


= lim 


bui — 


X x—foo 




x-s-oo bill 


= • • • = lim 


1-/ 






el- 

= 0 . 


Hence, e 


bpi xl-/3 o 

= C>(e"T^* ) = o{t-P). 


Lemma 3. Consider the following differential equation 

x{t) =— b{t))x{t), (11) 

where a(t) > 0, limt_»oo &(f) = 0 and K(A) > 0 (A £ C). For 
Ve > 0, there exists a positive constant M < 00 such that | 2 :(f)| < 
Me-(®(^)-")^‘o |a:(fo)|. 

Proof: The solution to lilt is x{t) = e ^‘^*0 
Let 5 = e/K(A) > 0. Then there must exist a finite positive constant 
T > to such that Vf > T, 1 — b{t) >1 — 5. 

• If to < r < L then 


|-Pi+i(t,to)| < [ a{r)\P.,^{r,to)\\Poit,r)\dr 

JtQ 

fa{r)dr. 

JtQ 

Let Mi+i = Mirjee, where e is the Euler’s number. Then it is 
calculated that 

|Pi+i(Lfo)| - 

^ Mi a{T)dr ^ A/p+ie^^-^4 

“ g(St(>')-7e)//„ “(slrfs g(9t(>‘)-F£)/t*(, 

Mi o(r)dr — Mi+ie ~^-^*0 

^(X(X)-^e) a(s)ds 

Define a function f(f) = Mi^ — Mi+ie~^. It is easy to see 
that maxj>o(/(0) = /(f) = 0. Therefore, |Pi+i(t,fo)| < 
Mi+ie~^^^^^ /to gy jgg mathematical induction, it 

can be proved that |Pi(f,fo)| < dto ^ j _ 

1, • ■ • , r, where M = max{l,Mi,-- - ,Mr-i}. Hence, there 
must exist a finite positive constant M such that ||$A(f,to)||i < 

Jf^a(s)ds ^ ^ 


< g-SS(>)//<*(s)(l-5)dOg-5R(>)/to <i(s)(l-'5+4-(>(s))do|^(j^)| 

• If to < f < r, then 


where M 2 = sup 4 g( 4 ^ j,. 


g®(^) /to a{3)(b{3)-S)ds 


< 00 . 


Let M = max{Mi,M 2 }. Then it is proved that x{t) 


< 


Lemma 4 (Lemma 2 in (Tj). Consider the following differential 
equation 

fA 1 -| 


where x{f) 
of ( 1121 is 


x(t) = —a(t) 


x(t), (12) 

r X r 


(a:i(f), ■ • ■ , Xr{t))'^ £ C. The state transition matrix 


Lemma 6. Consider the following differential equation 

x{f) = —a{t)Ax{t), (14) 

where a(t) > 0 and x{t) £ R". The solution to this differ¬ 
ential equation is x(f) = 'l?(f, fo)2;(fo), where is the 

state transition matrix. If all eigenvalues {Ai,--- , An} of A have 
positive real parts, then for Ve > 0, there exists a finite positive 
constant M such that ||'l>(f,fo )||2 < Me ‘>,(3)ds^ where 

■^min — min{SR(Ai)|i = 1, • • • , n} > 0. 

Proof: There exists a transformation matrix T such that 
P-Mr = A = diag(Ai, • • • , A,), where At £ (n £ N+ 

and ft = ti) is the Jordan block with the diagonal elements 

being At. 

The state transition matrix "(((fjfo) can therefore be written in the 
following form 

$(f, to) = T"Miag($Ai {t, to), ■■■ , $A, (I, to))T, (15) 

where thAt (^1 ^ 0 ) is defined by l ll3t . By Lemmaj^ there exists a finite 
positive constant Mi such that 

Therefore, ||$(f,fo )||2 < -^‘o where M — 

||r||2||T-i||2maxi{Mi}. ■ 

Lemma 7. Consider the following differential equation 


^\{t,to) 


'Po{t,to) P^Cto) 

0 PoiFto) 


0 0 


Pr-l{t,to)- 

Pr-2{t,to) 
Po{t,to) . 


x{t) = —a{t){A — B(t))x{t), (16) 

(■J 3 J where aft) > 0, x{t) £ R" and limt_»oo Bft) A a zero matrix. 
If all eigenvalues of A have positive real parts, then for Ve > 0, 
there exists a finite positive constant M such that ||^(f,fo )||2 < 
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M\e <^) !o o-Mda^ where is the state matrix of dl6t 

and Amin is defined in Lemma^ 

Proof: The solution to ([HI can be written as x{t) = 
$(t, to)a:(io) + a{s)^{t, s)B{s)x{s)d{s), where $(f,io) is the 

state transition matrix of 01. 

By Lemma 1^ there exists a finite positive constant Mi > 1 such 
that fo)|| < -^*0 Hence, 

\\x{t)\\2^<mt,to)h\\x{to)h 

+ f\{s)mt,s)h\\Bis)\\ 2 \\x{s)hdis) 

Jto 

Jto 

Consider another differential equation 


y{t) = -(Amin - e/2)a{t)y{t) + a{t)Mi\\B{t)\\ 2 y{t), (17) 


where y{t) € K. By Lemma there must exist a finite positive 
constant M 2 such that Mt > to, 

|t/(f)| < (18) 


— a (-^min e/2) a(s)ds 

Mi||i3(s)||2y(s)(is. Therefore, if 


The solution to JTTt is y(t) = e 

y{to) = ||a:(fo)||2, then for Vf > to, fo)a:(to)||2 = ||2;(f)||2 < 

Miy{t) < “^'’^‘*'’||a;(fo)||2- By the arbitrariness 

of x{to), there must exist a positive constant M < 00 such that 
Vf > to, ||T'(Lio)||2 < 


n-®) Sio °-Mda 


Lemma 8. Consider the following stochastic differential equation 

+ • • • + hii{t) + hoiit) = at), (19) 


where ^(f) is a mean square continuous random process. It is 
assumed that ^(f) is convergent to a random vector ((* in mean 
square, where -E||C*||2 < 00 . If all roots of polynomial s” + 
bn-is^~^ + ■ ■ ■ + bo = 0 have negative real parts, then we have 

(I) limt^oc £;||^(f) - C76o||i = 0 and limt^oc ||7‘H^)lli = 
0 , i = 1, • • • ,n. 

(II) IfEidt) - C) = 0(6"''*'") where a > 0 and fi & (0, 1), 

then E{at) - C/bo) = 0(e-'"‘7 and E{C"\t)) = 
C)(e-'^‘^), ,n. 

(III) If £;|C(i) - C*l^ = 0{t-^) where fi G (0,1), then 

E\ab) - = 0{t-^) and E\C^\t)\^ = OiM^), 

i = 1, ■ ■ ■ ,n. 

Proof: (I) see the proof of Lemma 2 in (H. 

(II) Let D denote the differential operator, namely D*^(f) = 
7*^(f)- Let {ri,-- - ,r„} denote the roots of polynomial s" + 
+ • • • + 60 = 0, where K(ri) < 0, i = I,-- - ,n. 
Then the stochastic differential equation 03 can be rewritten as 

nr=i(» - rm) = m- Let xi{t) = - rcm 

(/ = 1 , • • • , n - 1) and Xn{t) = ^(f). 

Part I; By the definition of a:i(f), it can be obtained that 
there must exist n — 2 constants fci,-- - , fc„_2 such that a:i(f) = 

2)"-7(t) + fcn-2D"-7(t) + ••• + kl^t) + nr=2(-^t)C(i). 

Therefore, according to (I), it can be obtained that xi{t) is conver¬ 
gent to n72(“’'*)C*/7 = —C*/ri in mean square. Furthermore, 
xi{t) — rixi{t) + <((/). Then, 


xi{t) = e'^^^^ ^°'^xi{to)+ ^\{C{s)-C) + C)ds, ( 20 ) 

Jto 


which follows that E{xi{f)} = Ri{t,to) + R 2 {t,to) + 
Rsityto), where Ri{t,to) = e’'i^‘“*“7i(fo), i?2(f,fo) = 

jt^ and Rs{t,to) = E{C)ds. 

It is easy to see that 

Ri{t,to) = 0{C^A- (21) 

Since i7(C(f) — C*) = ), there must exist a finite positive 

constant M such that |i7(^(f) — ^*)| < Therefore, 

|i?2(Lfo)| < M / ( 22 ) 

JtQ 

It is easy to see that limt_>oo = oo and 

ds = oo. Therefore, by L’HospitaTs rule, it is 

obtained that 
lim 

°° J to 

Hence jjjis together with ^T2\ 

leads to the fact that 

i? 2 (f,to) = 0 ( 6 "''*'’). ( 23 ) 


It can be calculated that 


RsAffo) 




0{e^A- 


(24) 


Therefore by ( 120b . ( 121b . ( 123b and ( 124b . it is proved that E{xi{t) — 
C*/(—ri)) = 0{e~^^ ). By the same procedure, it can be proved 
that 


E[xi{t)-C/l[]^,i-rj)) =0{e->^*''), ,n. (25) 

Part II: It is easy to see that ri7i(“’'0 = bo. Therefore, it is 
obtained from ( 125b that 


E{m - c/bo) = 0{e-^^ ). (26) 

By (|25j, i5(a:„_i(f) + r„((*/bo) = E{/{t)-fnCA + fnC/bo) = 
0(6“''*^), which together with ([26) leads to the fact that E^A)) = 

VnEiCt) - C/bo) + 0(e-'^^^) = 0(6-'“'’). 

Assume that there exists a positive integer k < n such 
that for Vi G {1,•••,/(:}, E{C'\t)) = 0(6“'"*'’). It is ob¬ 
tained from (|23 that E{xn-k-i{t) - n7n-fc(“’'t)C*/7) = 

E{UC^_A^-ri)at)-YY:=r.-A-r^)C/bo) = 0(6-'“"), which 
follows that EiC^+^At)) = 0(6-'^^'’) -f UC^-A-r^)E{m - 
C/bo) = 0(6-''* ). By the mathematical induction, it is proved 
that E{C'\t)) = 0(6-"*'’), i = ,n. 

(Ill) By (|^, it is obtained that 

C‘ " 

*1 H- < 3{R4,{t,to) + Rsityto) + Reitjto)), (27) 

ri 


where R4,{t,to) = |e’''^* *“7i(fo)|^, Rs{t,to) = 

|/;e’'*‘*"“’(C(7-C7dsf and Re{t,to) = \ + 

Cl|2 

ri I 

It is easy to see that 

EiRfft, to)) = 0(6^’*^’’*^*) = o(t-^). (28) 


According to the properties of mean square integral, it can be obtained 
that 

E{RCt,to))<(^j\^^^^^^^-^^E^\as)-Cfds'^ . (29) 
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Since E\({t) — there must exist two positive 

constants T and M such that for Mt > T, E\({t) — C,*\^ < Mt~^. 
Therefore, 


Jto Jto 


1 

X Ei{({s)-Cfds +je‘ 


*0 


s 2 ds. (30) 


k I ,2 \ 

nL E\C/h,-m\" 

j=i * " / 

By the mathematical induction, it is proved that (f) | = 

i = 1, • • • ,n. ■ 

References 


It is easy to see that 


JtQ 
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